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ABSTRACT 


In  connection  with  other  investigations  at  the  David  Taylor  Model  Basin , 
detailed  information  became  necessary  on  the  motion  of  air  bubbles  in  variable 
pressure  fields.  Since  no  information  on  the  subject  was  available,  a fundamen- 
tal study  of  the  motion  of  bubbles  was  undertaken.  As  an  initial  step,  experi- 
ments were  conducted  to  determine  the  drag  and  shape  of  sinj^e  air  bubbles 
rising  fireely  in  various  liquids. 

The  results  of  the  experiments  show  that  a complete  description  of  the 
motion  of  air  bubbles  is  not  possible  by  use  of  dimensionless  parameters  con- 
taining  the  usual  physical  properties  of  the  liquid  (viscosity,  surface  tension, 
density).  Three  types  of  bubble  shapes  were  observed  in  each  liquid,  namely 
spherical,  ellipsoidal,  and  spherical  cap.  For  a specific  liquid,  the  shape  of 
the  bubble  was  a functicm  of  its  volume. 

For  tiny  spherical  bubbles,  the  drag  coefficients  coincide  with  those  of 
corresponding  rigid  spheres.  With  increase  in  bubble  siae,  a decrease  in  the 
drag  as  compared  to  that  of  rigid  spheres  occurs  in  some  liquids.  Thus,  the 
drag  curves  of  the  spherical  bubbles  rising  in  various  liquids  fall  between  two 
limiting  curves,  namely  the  drag  curve  of  rigid  and  fluid  spheres,  respectively. 

It  was  not  possiUe  to  determine  a criterion  for  the  transition  of  the  bubbles 
from  **rigid"  to  fluid  spheres.  The  region  of  ellipsoidal  bubbles  extends  over 
different  ranges  of  Reynolds  numbers  for  the  various  liquids.  The  drag  coef- 
ficients of  spherical  cap  bubbles  are  independent  of  bubble  size  and  have  a 
constant  value  of  2.6. 

For  bubbles  (equivalent  radius  0.08  to  0.80  cm)  rising  in  tap  water  or 
in  water  containing  certain  surface-active  substances,  experiments  show  an 
increase  in  drag  as  compared  to  bubbles  in  pure  water. 

Results  of  tests  to  determine  the  effect  of  the  container  walls  on  the 
velocity  of  rise  are  presented.  A description  of  the  experimental  apparatus 
is  given.  A summary  of  the  theoretical  and  experimental  work  of  other  in- 
vestigators is  also  included. 

INTRODUCTION 

The  tests  described  in  this  report  are  a continuation  of  experiments  given  in  a previous 
Taylor  Model  Basin  report.  ^ These  experiments  were  initiated  in  conjunction  with  the  work 
under  projects  NS  718-201  and  NE  061-287.  The  present  tests,  continued  under  NS  715-102, 
were  ocmducted  for  the  purpose  of  investigating  the  motion  of  bubbles  rising  under  the 
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influence  of  gravity  as  an  initial  step  fbr  obtaining  information  on  the  behavior  of  bubbles  in 
variable  pressure  flelds. 

A body  rising  or  falling  under  the  influence  of  gravity  reaches  a constant  velocity 
(terminal  velocity)  when  all  forces  acting  on  it  are  in  equilibrium: 

Drag  * Buoyant  force  + Weight  - 0 

For  rigid  bodies,  the  drag  will,  in  general,  be  a complicated  function  of  the  geometry  of  the 
body,  the  velocity,  and  the  physical  properties  of  the  medium,  i.e.,  the  density  and  viscosity. 
For  fluid  bodies,  such  as  drops  and  gas  bubbles,  the  function  is  further  complicated  by  the 
fact  that  the  body  may  be  of  changeable  shape  and  that  properties  of  the  fluid  inside  the  glob- 
ule, such  as  density  and  viscosity,  and  interfacial  effects  may  also  be  important  factors.  In 
general,  the  shape  that  the  fluid  globule  assumes  is  some  complicated  function  of  the  hydro- 
dynamic,  viscous,  and  interfacial  forces. 

The  drag  of  fluid  bodies  may  either  be  equal  to  (as  is  the  case  for  small  bubbles)  or 
less  than  that  of  the  corresponding  rigid  body  depending  upon  the  conditions  at  the  interface. 
In  the  former  case  there  exists,  effectively,  a rigid  surface  at  the  interface;  in  the  latter  case, 
the  fluid  particles  at  the  boundary  have,  in  contrast  to  rigid  bodies,  nonvanishing  tangential 
velocities.  The  circulation  inside  the  fluid  body  thus  reduces  the  drag  of  the  body. 

The  experiments  described  in  this  report  consisted  of  the  determination  of  the  terminal 
velocity,  shape,  and  path  of  single  air  bubbles  rising  freely  in  various  liquids  as  a function 
of  bubble  size.  The  possible  effect  of  the  walls  of  the  container  on  the  velocity  of  rise  of 
the  bubble  was  also  investigated.  A summary  of  pertinent  theoretical  and  experimental  work 
of  other  investigators  is  included. 

THEORETICAL  SOLUTIONS 

Theoretical  solutions  for  the  drag  of  rigid  and  fluid  spheres,  moving  slowly  in  an  infi- 
nite medium,  have  been  obtained  for  the  following  boundary  conditions  at  the  surface  of  the 
sphere: 

1.  Rigid  spheres 

a.  Stokes’  solution^ 

(1)  Velocity 

(a)  The  velocity  vanishes. 

2.  Fluid  spheres 

a.  Hadamard-Rybczynski’s  solution**®’^ 

(1)  Velocity 

(a)  The  normal  velocity  component  vanishes. 

(b)  The  tangential  velocity  components  at  both  sides  of  the  surface  are  equal. 
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(8)  Strasses 

(a)  The  notmal  and  taogantial  atresseB  at  both  sides  of  the  surface  aie  equal.* 
b.  Boussinesq's  solution^ 

(1)  Velocity 


(a)  The  normal  velocity  component  vanishes. 

(b)  The  tangential  velocity  components  at  both  sides  of  the  surface  are  equal. 
(2)  Stresses 

(a)  The  normal  stress  at  the  inside  of  the  surface  is  larger  than  the  stress 
at  the  outside  due  to  the  dynamic  surface  tension.** 

(b)  The  tangential  stress  at  the  surface  is  increased  across  the  surface  due 
to  the  dynamic  increment  of  surface  tension. 

The  drag  of  a sphere  in  an  infinite  medium  of  uniform  velocity  V thus  becomes: 


1.  Rigid  sphere: 

2.  Fluid  sphere: 


D 

a.  D 


6wiifV 


2(1  3|t* 

3p  + 3^' 


b.  D 


6irpf  £/ 


« + r(2fi  f ifi') 
e + 3f(;t  + (i*) 


where  D is  the  drag, 

H is  the  coefflcioit  of  viscosity  of  the  medium, 
r is  the  radius  of  the  sphere, 

n'  ia  the  coefflcient  of  viscosity  of  the  fluid  inside  the  sphere,  and 
« is  the  coefflcient  of  surface  viscosity. 

Using  the  condition  of  equilibrium  for  a sphere  rising  under  the  influence  of  gravity, 
we  obtain  for  the  rigid  case: 

4 4 

6nnrV-  — nf^pS-  — wr^p'g 
3 3 

2 o 

^ “T  (p  - p')  Stokes*  Law 
^ #* 


•Tlw  prasaar*  incraaaa  acroaa  tha  aorfaca  dua  to  aurfaoa  tanaion  (>  waa  aaglactad  in  Hadamard'a  analjrain 
Ineluaton  of  thla  praaaiira  drop  in  tba  bonndary  condition  for  tha  normal  atraaa  doaa  not  chnnga  tiw  roaulta.  That 
ia  to  any.  anrfaca  tanaion  aa  manifaatad  only  in  a praaanra  incraaaa  inalda  tha  fluid  aphara  doaa  not  affact  ita 
motloa  (thia  raault  la  alao  obtalnad  by  putting,  in  Bouaainaaq'a  aaalyaia,  tha  coafflclant  of  aurfaca  viacoaity, 
aaa  anbaaqaant  footnota,  aquai  to  saro.) 

**Bouaalnaaq  aaaumad  that  a dynamic  aurfaca  tanaion  ealata  at  intarfacaa  in  motion.  Ita  magnituda  ia  givan  by 
tha  anm  of  tba  uaual  (atatie)  aurfaca  tanaion  and  tha  dynamic  ineramanb  Tha  dynamic  incramant  varlaa  ovar  th« 
aurfaca  of  tha  aphara  and  at  a givan  point  ia  proportional  to  tha  rata  of  dilatation  at  that  point.  Tha  conatant  of 
proportionality  ia  callod  (dim  to  ita  aimliarlt>  *o  tha  vlacoalty  coafflclant)  tha  coafllciant  of  aurfaca  viacoaity. 
(Surfaca  viacoaity  haa  tha  dlmanaiona  maaa/tima,  wblla  tha  dlmanalona  for  viacoaity  ara  maaa/langth  x tlma.) 
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where  p ie  the  density  of  the  medium  end  p' is  the  density  of  the  sphere. 

For  the  fluid  sphere: 

f/  - — ip  - p*)  .r*~  ^1*  Hsdamerd'Rybczynski’s  Lew 

® 3p'+2p 


V •*—  (p-  p*)  e + 2r  ill  II  ) Boussinesq's  Lew 

^ « + f(2p+  3p0 


For  the  cese  of  bubbles,  where  p"«  p end  p'«  p,  the  lest  two  expressions  reduce  to  t 


end 


I p 


V 


2 e + 3fp 
5"  p e + 2rp 


In  the  lest  equetion,  the  fector  epproeches  1 for  r epproeching  zero  or  e very 

€ + 2fp 

lerge;  it  epproeches  3/2  for  t « rp^  i.e.,  for  large  r or  for  e approaching  zero.  Hence, 
for  very  smell  bubbles,  Boussinesq's  solution  approaches  Stokes’  lew  as  a limit,  while  the 
other  limit  is  Hedamerd*Rybczynski's  solution. 

From  the  boundary  conditions  as  stated  above,  it  is  obvious  that,  for  the  Hedamard> 
Rybczynski  and  Boussinesq  solutions,  circulation  exists  inside  the  bubble.  For  the  Stokes 
solution,  of  course,  there  is  no  circulation  inside  the  sphere. 


DIMENSIONAL  ANALYSIS 

Since  an  analytical  solution  for  the  dreg  of  fluid  bodies  over  a large  range  of  sizes  is 
hardly  attainable,  dimensional  analysis  of  the  phenomenon  may  serve  to  correlate  the  experi* 
mental  results. 

In  the  case  of  fluid  bodies  the  following  physical  variables  are  usually  considered  as 
pertinent: 


V Velocity  of  the  body 

g Acceleration  due  to  gravity 

p Density  of  the  fluid  medium 

p'  Density  of  the  fluid  body 

I Length  parameter  of  the  fluid  body 
p Coefficient  of  dynamic  viscosity  of  the  fluid  medium 


s 


H*  Coeffloirat  of  dynamic  viaooaity  of  tbo  fluid  inaido  the  fluid  body 
9 Interfaoial  tension* 


The  complete  set  of  dimensionless  products  will  therefore  contain  five  such  products.  In 
principle,  it  is  immaterial  which  complete  set  is  chosen  for  the  representation  of  the  phenom- 
enon. For  example,  we  may  use 


Re,  We,  0 

M P 


f AC 0,  Re,  M,Re',^  -0 
P 


and  so  on,  whwe  Cp  is  the  drag  coefficient. 
Re  is  the  Reynolds  number 


We  is  the  Weber  number  (•^)- 
Re*  is  the  Reynolds  number  inside  the  fluid  body|« 

M is  a dimensionless  parameter 

\ pa’/ 

If  the  density  and  viscosity  of  the  gas  inside  a bubble  are  considered  negligible,  the  physical 
variables  are  reduced  to  six.  The  dimensionless  products  then  take  the  form  of 


and 


Re,  We)-0 
or 

Af)  -0 
or 

We,  M) 
etc. 


When  experimental  data  on  gas  bubbles  are  plotted  in  terms  of  dimensirmiess  products, 
complete  correlation  will  be  obtained  provided  the  variables  chosen  are  all  the  variables  upon 
which  the  phenomenon  depends. 

In  the  case  of  bubbles,  it  is  most  convenient  to  use  a length  parameter  which  is  based 
on  its  volume  rather  than  a physical  dimension  as  is  customary  for  rigid  bodies.  The  length 
parameter  ohosen  is  the  equivalent  radius  r,  where 


9/  volume 


•Tlw  eoeffleleat  of  aurfaM  viacoalty  la  nat  Inelndad  aiaoa  mara  la  ao  aaparlaianlal  avldaaoa  that  dmaale 
Mfaee  taaolaa,  as  peotslatsd  by  BaoBalsaaq;,  aslata. 
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For  bttbblM  rising  at  their  teminnl  velocity  the  dreg  coefficient  Cj^  can  then  be  written  as 

t'O- ; — 

instead  of  its  usual  form  of  f7/(l/2)p  A where  A is  the  projected  area.  For  the  analytical 
solutions  of  Stokes  and  HadamarckRybccynski  the  drag  coefficient  becomes,  respectively: 

C/,  “U  (Rigid  sphere) 

(Fluid  sphere) 

For  the  special  case  in  which  only  four  variaUes,  nsmely  the  velocity,  the  acoelera* 
tion  of  gravity,  the  density  of  the  fluid  medium,  and  the  equivalent  radius,  are  taken  as  per* 
tinent,  only  one  dimensionless  g^up,  the  <kag  coefficient,  is  obtained,  i.e., 

- constant 

This  solution  will  be  shown  to  apply  to  the  region  of  very  large  (spherical  cap)  bubbles. 

PREVIOUS  EXPERIMENTAL  WORK 
RATE  OF  RISE  OF  BUBBLES 

Interest  in  the  motion  of  air  bubbles  has  existed  for  many  years.  The  work  on  bubbles 
has,  however,  been  mostly  experimental  in  nature.  Exceptions  are  an  attempted  theoretical 
analysis  by  Theremin^  in  1829,  the  analytical  solution  of  Hadamard,^'^  Rybosynski,^  and 
Boussinesq,*'*  and  the  dimensional  analysis  of  Schmidt*  and  Rosenberg.* 

The  early  experimental  work  on  bubbles  was  largely  concerned  with  very  small  bubbles 
and  was  carried  out  for  the  purpose  of  determining  the  extent  of  Stokes  region.  Allen*  deter> 
mined  the  rate  of  rise  of  air  bubbles  in  water  and  in  aniline  up  to  bubble  radii  of  0.04  and 
0.06  cm,  respectively.  Arnold**  measured  velocities  of  small  air  bubbles  in  olive  oil  and  in 
aniline.  Bond  and  Newton**  investigated  eir  bobbles  in  syrup  and  in  water  glass  (sodium 
silicate). 

The  range  of  bubble  sixes  was  extended  by  other  investigators***  **  who  were  mainly 
interested  in  the  problem  in  ormnection  with  air*lift  pumps, *****  gas  absorption,**'**  or 
propagation  of  sound  in  liquids.**  These  expMiments  were  carried  out  in  water.  In  aubse* 
quent  years,  some  investigatimis  were  also  made  in  liquids  other  than  water.  Davies  and 
Taylor**  used  nitiobensmie  as  well  as  water  and  measured  velocities  of  large  bubblee. 
Temperley  and  Chambers**  extended  the  range  of  Taylor’s  experiments  in  water  to  bubUes 
of  equivalent  radii  up  to  approximately  6 cm.  Bryn**  made  tests  in  various  water>glyoerine 
and  water>ethyl  alcohol  mixtures.  Robinson**  measured  the  rate  of  rise  of  small  air  babbles 
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in  lubrionttng  oils.  His  results,  however,  show  oossiderable  scatter.  Reports  by  Piokert,^^ 
Pekeris,^^  Worster,^*  and  Datta  et  al  give  suminaries  of  the  results  of  experiments  of 
other  investigators. 

Qorodetskaya^^  investigated  the  effect  of  surfaoe-active  substances  on  the  rate  of 
rise  of  air  bubbles  in  water.  Further  tests  <»  air  bubbles  are  reported  in  References  81*86* 

In  addition,  a limited  number  of  tests  using  gas  bubbles  of  oxygen,  nitrogen,  and  a mixture 
of  carbon  dioxide  and  oxygen  have  been  carried  out  in  artificial  sea  water.^^'^*  A number  of 
tests  with  oxygen  bubbles  were  also  conducted  in  water^^and  in  aqueous  solutions  of  sodium 
hydroxide.^*  Recently,  Stuke^^  investigated  the  rate  of  rise  of  oxygen  bubbles  in  pure  (pro* 
snmably  distilled)  water  and  in  water  containing  surface*active  substances.  The  results  of 
the  tests  with  gas  bubbles  (given  in  the  Appendix)  show  no  significant  change  in  the  rate 
of  rise  of  the  bubbles  with  change  in  the  gas  inside  the  bubble. 

Because  of  the  scatter  of  previous  results  of  experiments  on  the  rate  of  rise  of  air 
bubbles  in  water,  Rosenberg^  repeated  these  tests  for  a large  range  of  air  bubble  sizes.  He 
showed  the  geometric  similarity  between  large  bubbles  of  spherical  cap  shape  and  suggested 
the  use  of  three  dimensionless  parameters,  the  drag  coefficient,  the  Reynolds  number,  and 
the  parameter  U for  describing  bubble  motion  in  liquids. 

WALL  EFFECT 

Previous  investigations  on  the  moticm  of  air  bubbles  in  liquids  were,  with  a few  ex* 
oepdons,  conducted  in  containers  of  limited  dimensions.  Only  for  very  large  bubbles  did 
Exner^^  and  Bryn^^  make  their  measurements  in  lakes.  Inasmuch  as  the  effect  of  the  walls 
of  the  container  on  the  rate  of  rise  of  bubbles  was  unknown,  it  was  generally  neglected. 
Miyagi^^  conducted  a few  tests  in  containers  of  different  sizes  and  found  that  a reduction  of 
4 percent  in  the  rate  of  rise  occurred  for  the  range  of  bubble  sizes  investigated.  Dubs^  ^ dc»- 
rived,  from  energy  considerations,  an  analytical  expression  for  the  wall  effect  and  concluded 
that  a bubble  of  the  same  radius  as  its  cylindrical  container  has  a velocity  of  rise  of  zero. 

It  is  clear  that  this  conclusion  is  in  error,  as  experiments  on  cylindrical  bubbles  have  shewn. 

As  indicated  previously,  no  analytical  solution  has  yet  been  obtained  for  flows  beyond 
the  region  of  slow  flow;  consequently,  the  much  more  difficult  problem  of  also  including  a 
finite  boundary  in  the  equations  of  motion  becomes  less  capable  of  solution.  For  very  slow 
flow  about  rigid  spheres  moving  in  an  infinite  cylindrical  c<»tainer,  Ladenburg^^  obtained 
an  analytical  solution  for  the  effect  of  the  boundary  on  the  drag  and,  crmsequently,  the  veloc- 
ity of  the  sphere.  McNown  et  al^^  arrived  experimentally  at  a wall  correction  coefficient  for 
rigid  spheres  descending  in  a cylindrical  container.  Since  air  bubbles  of  small  volume  rising 
in  water  behave  essentially  like  rigid  spheres,  this  correction  factor  may  be  applied  to  such 
bubbles  as  long  as  the  flow  is  still  in  the  Stokes  region. 

With  the  exception  of  a number  of  tests  for  large  bubbles,  no  data  concerning  the 
effect  of  the  boundary  on  the  rate  of  rise  of  gas  bubbles  beyond  the  region  of  slow  flow  are 
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available.*  However,  MOller^^  showed,  by  means  of  a dye  technique,  that  the  flow  about  a 
rigid  sphere  at  a lower  Reynolds  number  and  larger  boundary  dimensions  was  identical  to  that 
at  a higher  Reynolds  number  and  smaller  boundary  dimensions  and  therefore  that  the  effect 
of  the  walls  was  to  stabilise  the  flow  about  the  sphere.  These  results  for  rigid  spheres  be* 
yond  the  Stokes  region  of  flow  at  least  suggested  the  possibility  of  a similar  effect  for  the 
motion  of  gas  bubbles. 

CYUNDRICAL  BUBBLES 

At  this  point,  it  may  be  of  interest  to  mention  experiments  on  a special  form  of  finite 
boundary  dimensions,  i.e.,  the  case  of  cylindrical  bubbles.  This  term  was  first  used  by 
Qibson^^  and  applies  to  the  type  of  bubble  formed  when  a long  cylindrical  tube  filled  with 
liquid  is  emptied  from  below  or  when  a large  amount  of  air  is  introduced  through  the  bottom 
of  the  tube. 

Qibson  investigated  the  velocity  and  shape  of  these  bubbles  in  water.  Ward  and  Kess* 
ler^^  conducted  tests  in  pipes  of  various  diameters.  Hattori^^  was  interested  in  the  problem 
in  connection  with  the  possibility  of  evaluating  the  surface  tension  of  a liquid.  Hence,  he 
was  concerned  with  tubes  of  small  diameter,  since  the  so-called  critical  tube  diameter  (below 
which  the  bubble  no  longer  rises  but  remains  stationary)  is  a function  of  the  surface  tension 
of  the  liquid.  Dumitrescu^^  obtained  an  analytical  expression  for  the  velocity  of  a cylindrical 
bubble  by  neglecting  viscous  and  surface  tension  forces,  thus  reducing  the  problem  to  one  of 
potential  flow.  The  differential  equation  for  the  velocity  potential  together  with  the  existing 
boundary  condition  yields  a solution  for  the  velocity  of  rise  as  a function  of  the  tube  diameter 
only.  His  experimental  tests  in  water  show  that  for  a tube  of  sufficiently  large  diameter 
(3  cm  for  water  at  room  temperature)  the  measured  velocities  agree  very  closely  with  the 
theoretical  values.  Therefore,  for  large  bubbles,  the  physical  properties  of  the  liquid  no 
longer  have  any  effect  on  the  flow  about  the  bubble  and  the  bubbles  are  geometrically  similar. 
Davies  and  Taylor^^  investigated  the  shape  and  rate  of  rise  of  cylindrical  bubbles  in  order 
to  obtain  a better  understanding  of  the  pressure  distribution  of  spherical  cap  bubbles  in  an 
infinite  medium. 


SCOPE  OF  THE  PRESENT  INVESTIGATION 

The  present  investigation  was  initiated  in  connection  with  a program  of  study  of  the 
behavior  of  air  bubbles  in  water  at  variable  pressure  gradients.  Since  extensive  experimenta- 
tion was  required  for  direct  experimental  study  of  the  motion  of  air  bubbles  in  such  pressure 


*A  pw>  by  Coppock  and  Maiklejobn^  haa  racantly  coma  to  the  attantion  of  tha  autbora.  From  taata  con* 
ductad  with  air  bubblaa  in  watar,  tbay  concluda  that  no  wall  affact  axiata  for  bubblaa  ranging  in  aqulvalant 
radiua  from  0.01  to  0. 1 cm  rlaing  In  a tuba  of  5 cm  dlamatar. 
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fields,*  an  sltemste  approach  appeared  more  feasible.  It  consists  of  calOulating  the  motion 
of  the  bubbles  in  variable  pressure  fields  from  a knowledge  of  bubble  drag  at  various  constant 
pressure  gradients.  Experimental  data  on  the  drag  of  air  bubbles  at  various  pressure  gradients 
are  essential  in  this  procedure.  However,  only  data  on  the  motion  of  bubbles  in  pressure  gra* 
dients  produced  by  gravity  were  available.  Therefore,  information  on  bubble  motion  in  water 
at  pressure  gradients  other  than  gravity  became  necessary.  This  information  could  be  obtained 
by  investigating  the  rise  of  bubbles  in  liquids  having  the  same  physical  properties  as  water 
with  the  exception  of  the  density.  With  all  other  properties  of  the  liquids  identical,  varying 
the  density  would  be  equivalent  to  varying  the  pressure  gradient.  This  approach,  however,  is 
not  practicable,  since  there  are  no  liquids  available  which  possess  such  properties.  The 
other  approach  is  to  investigate  the  rise  of  air  bubbles  in  various  liquids  having  different 
physical  properties  and  then  to  attempt  to  correlate  the  results  in  terms  of  nondimensional 
parameters.  The  available  information  on  the  rise  of  bubbles  in  different  liquids  was  too 
meager  to  allow  definite  conclusions  regarding  the  significance  of  the  parameters  suggested 
by  Reference  1.  The  present  investigation  was  therefore  initiated  with  the  purpose  of  deter* 
mining  the  nondimensional  parameters  for  bubble  rise  by  investigating  bubble  motion  in  a 
number  of  liquids  of  different  physical  properties.  If  it  were  found  that  the  motion  of  air  bub- 
bles rising  freely  in  a liquid,  that  is  to  say  the  motion  in  the  pressure  field  produced  by  gravi- 
ty, could  be  described,  for  example,  in  terms  of  the  drag  coefficient,  the  Reynolds  number, 
and  the  parameter  i/(-  the  results  thus  obtained  could  be  used  in  evaluating  the 

drag  of  bubbles  in  water  at  pressure  gradients  other  than  that  produced  by  gravity.  To  do 
this  it  would  have  to  be  shown  that  the  nondimensional  parameters  used  for  the  freely  rising 
bubbles  are  also  applicable  to  other  pressure  fields.  This  might  be  accomplished,  for  example, 
by  comparing  the  results  of  a bubble  experiment  in  water  in  a nongravity  pressure  gradient 
with  the  results  of  bubbles  rising  freely  in  various  liquids  at  identical  "/V”  number. 

By  conducting  the  tests  on  the  rise  of  bubbles  in  various  liquids  in  a large  tank,  the 
possibility  of  the  effect  of  the  tank  walls  on  tiie  velocity  of  the  bubbles  is  eliminated.  Since 
the  high  cost  of  many  desirable  liquids  makes  the  use  of  a large  tank  impractical,  it  became 
necessary  to  determine  the  possible  effect  of  the  walls  on  the  velocity  over  the  range  of  bub- 
ble sixes  to  be  tested.  This  investigation  consequently  acquired  two  purposes: 

1.  The  determination  of  the  effect  of  variation  of  liquid  properties  on  the  motion  of  air 
bubbles. 

2.  The  evaluation  of  wall  effect. 


*Ezploratory  •sparioMnts  of  such  • utnro  are  rapoitod  in  Roforenc*  35. 

••Thin  pwamntor  la  glvan  In  a mora  ■anaral  form  aa  {p*  Pp)/(p*<r®)  (whara  Vp  la  the  praaanra  grndlant;  for  a 
gravity  flald  Vp  "•  Tharafora,  for  a apaclflc  liquid,  it  la  propoitlonai  to  tha  praaanra  gradlant 
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EXPERIMENTATION 

The  experimenUl  study  consisted  of  messuring  the  terminsl  velocity  of  individual  bub* 
bles  of  various  sises  rising  in  eight  liquids.  It  also  included  the  determination  of  the  effect 
of  the  walls  of  the  container  on  the  bubble  velocity.  The  experimental  apparatus  consisted 
of  test  tanks,  means  for  determining  the  physical  properties  of  the  liquids,  a regulated  bubble 
supply,  and  means  for  measuring  bubble  size  and  velocity.  Details  of  the  experimental  appa- 
ratus and  procedure,  the  generation  of  the  bubbles,  and  the  test  liquids  are  given  in  the  fol- 
lowing paragraphs. 

TEST  TANKS 

In  order  to  obtain  free  bubble  rise  and  to  reduce  the  effects  of  such  boundaries  as  the 
bottom  of  the  tank  and  the  free  liquid  surface,  the  containers  for  the  liquids  had  to  be  at 
least  2 ft  high.  The  tests  were  performed  in  three  transparent  wall  tanks;  the  large  one  was 
of  3 X 8 ft  cross  section  and  5 ft  height,  the  medium  one  was  of  1 x 1 ft  cross  section  and 
3 1/2  ft  height,  and  the  small  one  was  of  6 x 6 in.  cross  section  and  2 ft  hei^t.  In  addition, 
tests  were  also  performed  in  an  insert  of  6 x 6 in.  cross  section  and  20  in.  height,  placed  in 
the  center  of  the  medium  tank. 

Since  the  large  tank  was  of  sufficiently  large  dimensions,  no  significant  wall  effect 
was  expected.*  The  dimensions  of  the  medium  tank  were  chosen  large  enough  to  reduce  wall 
effects,  yet  small  enough  to  allow  use  of  a variety  of  liquids.  The  small  tank  and  the  insert 
provided  an  additional  tank  size.  It  was  intended  that  if  wall  effect  existed,  the  results  ob- 
tained in  the  finite  containers  would  be  extrapolated  to  the  case  of  an  infinite  medium. 

TEST  LIQUIDS 

The  eight  test  liquids  used  were  water  (at  three  different  temperatures),  Varsol,** 
methyl  alcohol,  turpentine,  water  containing  0.42  percent  (by  volume)  Glim,***  mineral  oil, 
and  two  com  syrup-water  mixtures.  Turpentine  was  selected  as  one  of  the  test  liquids  be- 
cause at  room  temperature  it  has  the  same  viscosity  as  cold  water.  One  of  the  com  syrup 
mixtures  had  approximately  the  same  viscosity  as  the  mineral  oil  (see  Table  1). 

The  viscosity  of  the  liquids  was  measured  by  means  of  ordinary  and  modified  Ostwald 
viscosimeters.  The  accuracy  of  measurement  of  viscosity  was  1.5  percent  and  0.5  percent. 


*A  few  of  the  previoue  Taylor  Model  Baain  teata  ^ arere  repeated  in  the  large  tank  to  obaarve  any  change  In 
reaulta.  Theae  prealoua  teata  ware  conducted  in  a tank  of  4 1/2  X 25  ft  croea  aectlon  and  9 ft  height  with  8 ft 
depth  of  filtered  water  at  room  temperature,  naing  one  end  of  the  tank  for  the  teata. 

**A  trade  name  (Standard  Oil  Company)  for  mineral  apirlta  ptaavy  naphtha),  a petroleum  dlatiUata. 

***Gllm  (Antarox  A-480),  a aurfaca-actire  agent,  la  the  trade  name  (3.T.  Babbitt,  Inc.)  of  a non-ionic,  liquid 
detergent,  a condensation  product  of  ethylene  oxide  and  lauryl  alcohol. 
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TABLE  1 


Summary  of  Liquid  Properties 


Surface 

"ir 

Temperature 

Viscosity 

Density 

Tension 

Liquid 

deg  C 

#* 

poises 

P 

gm/cc 

<r 

dynes/cm 

ilL 

pa^ 

Wat«r 

19 

0.0102 

0.998 

72.9 

0.26  X ir^o 

Water 

21 

0.0098 

0.998 

72.6 

0.24  X 10"^® 

Cold  Water 

6 

0.0147 

0.999 

74.8 

1.08  X ir*® 

Hot  Water 

49 

0.0056 

0.989 

68.1 

0.307  X 

Glin  Soiution 

19 

0.0103 

1.000 

32.8 

2.78  X 10"“ 

Mineral  Oil 

27.5 

0.580 

0.866 

20.7 

1.45  xir2 

Varsol 

28 

0.0085 

0.782 

24.5 

4.3  xir“ 

Turpentine 

23 

0.0146 

0.864 

27.8 

24.1  xir“ 

Methyl  Alcohol 

30 

0.0052 

0.782 

21.8 

0.89  X ir“ 

62  percent  Corn  Syrup 
and  Water 

22 

0.550 

1.262 

79.2 

0.155  X ir^ 

68  percent  Corn  Syrup 
and  Water 

21 

1.090 

1.288 

79.9 

0.212 X ir2 

56  percent  Glycerine 
and  Water  (Bryn^^) 

18 

0.0915 

1.143 

69.9 

1.75  X ir^ 

42  percent  Glycerine 
and  Water  (Bryn^^) 

18 

0.043 

1.105 

71.1 

4.18  X 10t8 

13  percent  Ethyl  Alcohol 
and  Water  (Bryn^) 

22 

0.0176 

0.977 

43.5 

1.17  X lor* 

Olive  Oil  (Arnold^®) 

22 

0.73 

0.925 

34.7 

0.716  X ir^ 

Syrup  (Bond'^) 

17 

180 

1.48 

91 

0.92  X 10« 

respectively.  The  surface  tension  was  determined  by  the  capillary-rise  method  (accuracy  of 
measurement:  3 percent)  and  the  sped  He  gravity  of  the  liquids  was  obtained  by  means  of 
hydrometers  (accuracy  of  measurement:  1 percent).  These  physical  properties  were  measured 
following  the  completion  of  each  test.  They  are  summarised  in  Table  1 together  with  those  of 
liquids  used  by  several  other  investigators. 

For  Varsol  and  water  (room  temperature  and  hot),  tests  were  conducted  in  all  three 
tanks;  for  cold  water  and  mineral  oil,  in  the  medium  tank  and  insert;  and  for  all  other  liquids, 
in  the  small  tank  only. 
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BUBBLE  GENERATION 

Small  bubbles  were  generated  by  means  of  hypodermic  needles  and  glass  nozzles  of 
various  sizes.  The  larger  bubbles  were  obtained  by  use  of  a dumping  cup,  which  was  inverted 
to  release  the  air  bubble.  The  nozzles  and  needles  were  connected  to  a brass  tube  which 
was  fastened  to  a sliding  mechanism  (Figure  1).  This  sliding  mechanism  allowed  the  tips  of 
the  various  nozzles  to  be  placed  at  the  identical  position.*  The  air  was  supplied  from  a com- 
pressed air  bottle.  A needle  valve  regulated  the  air  flow  so  that  bubbles  were  released  at  the 
interval  desired. 

DETERMINATION  OF  BUBBLE  SIZE 

The  bubble  size  was  determined  by  "weighing”  a sufficient  number  of  bubbles  in  the 
inverted  funnel  (Figure  1)  by  means  of  an  analytical  balance.  Since  the  density  of  air  is 
negligible  in  comparison  to  that  of  a liquid,  the  difference  in  balance  reading  equals  the 


•Thla  davic*  •llmlMtad  focuBinc  of  tha  camora  altar  aanh  changa  of  noaala.  Hia  caoiara  waa  uaad  to  datar> 
adaa  volocitjr,  path,  and  ahapo  of  tba  bubbla. 
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buoyancy  of  the  bubbles  (i.e.,  it  equals  the  volume  of  the  bubbles  times. the  density  of  the 
liquid).  The  change  in  balance  reading  was  always  at  least  0.2  gm,  resulting  in  an  accuracy 
of  measurement  of  volume  of  1 percent.  The  volume  of  the  individual  bubble  was  obtained  by 
dividing  the  total  volume  by  the  number  of  bubbles  collected  in  the  funnel.  A comparison  of 
photographs  of  different  bubbles  showed  that  the  bubble  size  did  not  vary  if  the  frequency  of 
bubble  generation  remained  constant.  Large  bubbles  from  the  dumping  cup  were  weighed 
individually.  The  volume  of  the  bubble  was  adjusted  fbr  the  change  in  pressure  due  to  differ* 
ence  in  depth  between  the  level  at  which  the  rate  of  rise  is  determined  and  the  level  of  the 
inverted  funnel.  This  was  done  by  use  of  the  general  gas  law  at  constant  temperature,  taking 
into  account  the  partial  pressure  of  the  saturated  vapor  at  test  temperature.* 

Tiny  spherical  bubbles  could  not  be  generated  at  a frequency  to  allow  a sufficient 
number  to  be  collected  in  the  funnel,  hence  their  size  was  determined  from  the  photographic 
record.  No  correction  for  change  in  depth  is  then  needed. 

OTHER  EXPERIMENTAL  TECHNIQUES 

To  avoid  any  changes  in  the  volume  of  the  bubble  due  to  air  interchange  with  the 
liquid,  the  latter  was  saturated  with  air  prior  to  actual  testing.  This  was  accomplished  by 
stirring  the  liquid  and  by  blowing  air  through  it. 

The  temperature  of  the  liquids,  with  the  exception  of  cold  and  hot  water,  was  room 
temperature,  which  varied  little  throughout  the  day.  Water  was  cooled  by  circulation  through 
a water  cooler;  it  was  heated  with  immersion  heaters  or  obtained  directly  from  the  hot  water 
faucet.  Both  Rltered  and  tap  water  were  used  in  the  tests. 

Uniformity  of  liquid  temperature  was  achieved  by  means  of  mechanical  stirring  before 
each  test.  Frequent  checks  of  temperature  at  various  locations  inside  the  tank  were  made  by 
means  of  immersion  thermometers.  In  the  process  of  stirring,  small  bubbles  appeared  in  the 
liquid.  The  irregularity  of  motion  of  these  small  bubbles,  which  were  still  present  after  com* 
pletion  of  the  stirring,  served  as  an  indication  of  the  presence  of  residual  turbulence  in  the 
liquid.  In  sufficient  time,  the  motion  of  the  small  bubbles  always  became  regular  and  hence 
indicated  that  the  residual  turbulence,  if  sUll  present,  was  not  large  enough  to  affect  the 


*D«taila  of  thia  corroctlon  ara  aa  followa; 


whara  la  tha  aoluma  of  tha  bubblo  at  tba  canara  laval, 

Vq  la  tha  Toluaia  of  tha  bubbla  aa  datarminad  by  walgbing, 
la  tba  abaoluta  praaaura  at  tha  funaal, 
la  tba  abaoluta  praaaura  at  tha  camara  laval,  and 
p'  la  tba  vapor  praaaura  of  tha  liquid  at  tha  taut  tamparatura. 
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motion  of  the  bubbles.  The  actual  tests  were  not  begun  until  all  the  small  bubbles  reached 
the  surface  of  the  liquid. 

The  rate  of  bubble  flow  was  then  regulated  by  the  needle  valve  so  as  to  release  bubbles 
with  a minimum  spacing  of  24  in.  This  reduced  tlie  effect  of  the  wake  created  by  the  passage 
of  a bubble  on  the  motion  of  a bubble  following.*  At  higher  bubble  rates,  tiie  velocity  of  the 
individual  bubble  is  increased.  

The  same  precaution  was  observed  for  the  larger  bubbles  that  were  formed  by  dumping. 
An  additional  precaution  was  to  rotate  tlie  dumping  cup  with  steady  speed  in  order  to  avoid 
splitting  of  the  bubble  or  the  formation  of  satellites  upon  release.  The  slow  passage  of  the 
air  through  the  brass  tube  inside  the  tank  allowed  the  air  to  reach  tlie  temperature  of  the  li- 
quid. Contact  of  the  air  at  the  notale  tip  or  inside  the  dumping  cup  with  the  liquid  allowed 
saturation  of  the  air  with  liquid  vapor,  so  that  the  air  bubble  can,  in  each  instance,  be  as- 
sumed to  be  saturated  with  the  vapor  of  the  liquid  in  which  it  rises. 

MOTION  PICTURES  AND  THEIR  EVALUATION 

The  velocity,  path,  and  shape  of  the  bubbles  were  obtained  from  motion  pictures  made 
with  a Mitchell  86mm  camera  using  a special  lens  attachment  to  permit  close-ups.  Film  speeds 
of  25  to  35  frames  per  second  and  back  lighting  horn  a white  reflector  were  used.  For  the  first 
few  tests,  the  film  speed  was  obtained  by  photographing  a rotating  clock  dial;  subsequently  a 
neon  timing  light  with  a 60-cycle  voltage  source  was  utilized.  The  Him  speed  was  determined 
from  the  marks  of  the  timing  light  on  the  film.  The  field  of  the  camera  varied  from  1.4  x 1.8 
to  1.75  X 2.3  in.  depending  upon  the  refractive  index  and  horizontal  depth  of  liquid.  A trans- 
parent scale  photographed  in  the  plane  of  tlie  bubble  provided  the  distance  scale  factor  for 
the  evaluation  of  displacement  and  size.  The  camera  lens  was  placed  at  approximately  the 
midpoint  between  the  liquid  level  and  the  bottom  of  the  tank  for  all  tests.  The  camera  loca- 
tion was  in  each  instance  sufficiently  above  the  nozzle  tip  so  that  the  bubbles  reached  their 
terminal  velocity  before  passing  in  front  of  the  camera.  A summary  of  camera  location  and 
depth  of  liquids  in  the  tanks  is  given  in  Figure  2. 

Changes  in  bubble  volume  due  to  differences  in  liquid  depth  were  minimized  by  making 
velocity  measurements  over  a very  short  vertical  displacement  (less  than  2 1/2  in.).  The 
rate  of  rise  of  bubbles  was  determined  by  measuring  the  displacement  of  a bubble  from  a 
reference  point  on  successive  frames  of  the  Him  by  means  of  a Bausch  and  Lomb  contour- 
measuring  projector  using  a magnification  of  twenty-five  (Figure  3).  These  displacements 
were  then  plotted  against  the  frame  number.  The  straight-line  plot  indicates  that  the  veloc- 
ity of  the  bubble  remained  constant  during  the  time  it  passed  the  field  of  Uie  camera.  From 
the  slope  of  the  line,  the  frame  speed,  and  the  scale  factor,  tiie  velocity  of  the  bubble  is 
computed. 

«Kapi*r^  Bhowcd  th«  abaanca  of  proximity  affact  for  air  bubblaa  in  water,  ranging  in  equivaiant  radiua  from 
0. 14  to  0,38  cm,  if  the  frequency  waa  balow  30  bubblaa  par  minute. 


2<rl-  4- ■ 3.33 


Liquid 

Av.  Temp, 
deg  C 

Tank 

Width  of  Tank 

W 

inches 

1 

Height  of  Liquid 
h 

inches 

% 

Height  of  Liquid 
above  Camera 
hi 

inches 

Water  (Fiitered) 

19 

Large 

36 

48 

1.33 

22.5 

Water  (Tap) 

20 

Large 

36 

48 

1.33 

25 

22 

Medium 

12 

31.5 

2.62 

15.5 

20 

Insert 

6 

36.5 

3.33 

20 

Water  (Filtered) 

6 

Large 

36 

48 

1.33 

23 

5 

Medium 

12 

30 

2.50 

14 

Water  (Tap) 

48 

Large 

36 

49 

1.36 

23.5 

49 

Medium 

12 

39 

3.25 

22.75 

48 

Small 

6 

19 

3.17 

9 

Gliai  Solution 

19 

Small 

6 

20 

3.33 

9.5 

Mineral  Oil 

28 

Medium 

12 

36 

3.00 

19 

27 

Insert 

6 

36.5 

3.33 

19.5 

Varsol 

28 

Large 

36 

39.5 

1.10 

24.5 

28 

Medium 

12 

38.5 

3.21 

21.5 

25 

Insert 

6 

36 

3.33 

19 

Turpentine 

23 

Small 
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Figure  2 - Test  Conditions 
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RESULTS 

TERMINAL  VELOCITY  OF  AIR  BUBBLES 

Th«  results  of  tests  to  detemine  the  velocity  of  rise  of  air  bubbles  in  various  liquids 
are  most  oonvenimtly  presented  as  a function  of  the  equivalent  radius  of  the  bubble,  defined 
as  the  radius  of  a sphere  having  the  same  volume  as  the  bubble.  Figures  4«18  show  the  ter- 
minal velocity  of  air  bubbles  rising  freely  in  tap  (unfiltered)  and  in  filtered  water  (including 
data  ftom  other  investigators),  in  water  containing  Glim,  in  mineral  oil,  Varsol,  turpentine, 
methyl  alcohol,  and  two  com  syrup-water  mixtures  as  a function  of  the  equivalent  radius. 
Figure  14  presents  Bryn's  results  in  an  ethyl  alcohol-water  mixture  and  two  glycerine-water 
mixtures.*  Figure  15  summarizes  all  velocity  curves  (except  those  for  tap  water).  A compila- 
tion of  the  properties  of  the  liquids  is  given  in  Table  1 (see  page  11). 

In  general,  the  results  as  seen  from  Figure  15  indicate  that  for  small  (spherical)  air 
bubbles  of  given  volume,  the  viscosity  of  the  liquid  is  the  most  important  property  determining 
the  rate  of  rise.  Very  large  bubbles  (spherical  caps)  rise  independently  of  the  properties  of 
the  liquid. 

WALL  AND  PROXIMITY  EFFECTS 

As  indicated  previously,  the  effect  of  the  container  walls  on  the  velocity  of  a bubble 
had  to  be  determined  if  the  reeults  of  tests  conducted  in  a tank  of  limited  dimensions  were 
to  be  applied  to  bubble  motion  in  an  infinite  medium.  Tests  were,  therefore,  conducted  in 
tanks  of  different  sizes  in  water,  Varsol,  and  mineral  oil.  Figure  4 gives  the  results  of  tests 
conducted  with  filtered  water  by  several  investigators  including  the  Taylor  Model  Basin. 

The  cross  sections  of  the  containers  used  are  also  indicated  in  the  figure.  No  wall  effect 
is  noticeable  from  these  results.  For  example,  Gorodetskaya's  results  for  bubbles  ranging 
from  0.01  to  0.07  cm  rising  in  a tube  of  5 cm  diameter  show  no  wall  effect  when  compared 
with  results  of  tests  conducted  in  larger  containers.  The  results  of  the  present  experiments, 
given  in  Figures  5,  6,  7,  9,  and  10,  show  within  experimental  accuracy,  the  absence  of  any 
wall  effect  tor  the  range  of  bubble  sizes  tested.  Subsequent  tests  in  the  other  liquids  were 
made  in  the  small  tank  only  and  the  results  may  be  applied  to  the  case  of  an  infinite  medium. 

No  systematic  investigation  was  made  of  vertical  proximity  effect,  i.e.,  the  effect  of 
the  wake  created  by  the  passage  of  a bubble  on  the  motion  of  a bubble  rising  at  a distance 


MaalUi  of  tbo  81  poreoat  (by  wolght)  glycorlM-wator  mixturo  havo  boon  omittod.  Biyn  prooontod  thooo 
roonltn  in  toraw  of  drof  eeofficloat  and  Roynoldo  nunobor,  from  which  tho  tormlnal  roloclty  can  bo  computod. 

In  tho  rogton  of  bobbla  also  wharo  tho  rata  of  rlao  la  shown  to  bo  a fOnctlon  of  also  only  (banco  a common  voloc- 
Ity  eofvo  for  all  Uqulda;  ooo  FIzom  15).  tho  volodty  curro  fbr  tho  81  porcont  mlsturo  falla  appMclably  abovo 
tho  ooo»on  enrva.  Tho  dlocroponey  to  probably  duo  to  orronoona  oralnatlon  ot  tho  two  dlmonolonlooo  paramotara 
for  mo  81  poreont  glyoortno-wntor  mtattofo. 
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Figure  4 - Terminal  Velocity  of  Air  Bubbles  in  Filtered  or  Distilled 
Water  as  a Function  of  Bubble  Size 


Figure  6 - Terminal  Velocity  of  Air  Bubbles  in  Tap  Water 
as  a Function  of  Bubble  Size 
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Figure  6 • Terminal  Velocity  of  Air  Bubbles  in  Cold  Filtered  Water 
as  a Function  of  Bubble  Site 
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Figure  7 - Terminal  Velocity  of  Air  Bubbles  in  Hot  Tap  Water 
as  a Function  of  Bubble  Site 
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Figure  10  • Terminal  Velocity  of  Air  Bubbles  in  Varsol 
as  a Function  of  Bubble  Site 


Figure  11  - Terminal  Velocity  of  Air  Bubbles  in  Turpentine 
as  a Function  of  Bubble  Size 
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Figure  18  • Terminal  Velocity  of  Air  Bubbles  in  Com  Syrup*Water  Mixtures 
as  a Function  of  Bubble  Sice 
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Figure  14  * Terminal  Velocity  of  Air  Bubbles  in  Liquids  as  Obtained 

from  Bryn's  Data 


Figure  16  - Comparison  of  Terminal  Velocities  of  Air  Bubbles  in  Various  Liquids 
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below.  Suoh  effects  were  avoided  in  the  experimoits  by  sufficient  spacing  between  the  bub* 
bles.  However,  the  results  of  a few  special  observations  indicate  that  proximity  effects  may 
be  appreciable.  For  example,  tests  in  mineral  oil  show  an  increase  of  9 percent  and  89  per* 
cent  for  bubbles  of  equivalent  radius  of  0.17  cm,  rising  7.7  cm  and  8.2  cm  apart,  respectively. 
Napier^®  observed  an  increase  of  8 percent  for  bubbles  of  0.14  cm  radius  in  water,  rising 
6 cm  apart.  The  presence  of  the'  wake  in  the  liquid  thus  results  in  higher  velocities  of  rise 
of  the  bubble. 

NONDIMENSIONAL  PRESENTATION  OF  BUBBLE  DATA 

In  a previous  section  it  was  pointed  out  that  presentation  of  tiie  experimental  data  on 
air  bubbles  in  terms  of  dimensionless  products  gives  complete  correlation  provided  the  vari- 
ables considered  in  the  analysis  are  complete  and  pertinent.  The  results  of  the  Taylor  Model 
Basin  bubble  tests  and  those  of  Arnold,^®  Bond  and  Newton,^*  and  Bryn^^  are  given  in  terms 
of  the  drag  coefficient,  Reynolds  number  and  the  parameter  M in  Figure  16.*  Figure  17  pre* 
sents  the  bubble  data  in  terms  of  the  drag  coefficient,  Weber  number,  and  the  parameter  M. 

The  curve  for  filtered  or  distilled  water  at  a temperature  of  19  deg  C was  drawn  through  points 
obtained  from  the  experiments  of  Bryn  and  the  Taylor  Model  Basin  tests. 

Examination  of  Figure  16  or  17  shows  no  systematic  arrangement  of  the  curves  with 
change  in  the  parameter  A/,  wtiich  is  constant  for  a specific  liquid.  It  can,  therefore,  be  con- 
cluded that  neither  of  the  nondimensional  sets  presented  nor  any  other  complete  set  using 
the  same  six  variables  (namely,  velocity,  acceleration  of  gravity,  density  and  viscosity  of 
the  liquid,  surface  tension,  and  equivalent  radius)  is  sufficient  for  a complete  description  of 
bubble  motion. 

The  question  now  arises  whetiier  correlation  of  bubble  data  could  be  obtained  by  using 
two  additional  dimensionless  parameters,  for  example,  the  liquid  to  air  viscosity  and  density 
ratios  or  the  Reynolds  number  inside  the  bubble  and  the  density  ratio,  etc.  The  results  of 
the  experiments  conducted  do  not  permit  conclusions  regarding  the  importance  of  these 
parameters.  A short  discussion  of  the  significance  of  the  internal  Reynolds  number  will  be 
given  in  subsequent  sections. 

SPHERICAL  BUBBLES 

It  was  observed  in  the  experiments  that,  as  the  bubble  size  was  increased,  a change 
in  bubble  shape  from  spherical  to  ellipsoidal  to  spherical  cap  shape  occurred  in  all  liquids. 
Very  small  bubbles  are  spherical.  Larger  bubbles  are  flattened,  i.e.,  ellipsoidal  in  shape, 
whereas  very  large  bubbles  assume  a spherical  cap  shape.  Of  course,  the  volumes  at  which 


*Th*  results  for  tsp  wstsr  snd  for  wstsr  contsinlna  Glim  sie  not  shown.  They  will  be  discussed  in  subsequent 
sections. 
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Figure  16  - Drag  Ck>efficient  as  a Function  of  Reynolds  Number  for  Air  Bubbles 
Rising  at  Their  Terminal  Velocity  in  Various  Liquids 
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Figure  17  - Drag  Coefficient  as  a Function  of  Weber  Number  for  Air  Bubbles 
Rising  at  Their  Terminal  Velocity  in  Various  Liquids 
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tli«8e  transitions  occur  vary  with  the  liquids.  Photographs  of  typical  shapes  are  shown  in 
Figure  18.  It  should  be  noted  here  that  some  of  the  shapes  shown  in  these  photographs  are 
instantaneous  shapes,  since  the  shape  of  large  bubbles  does  not  remain  constant  during  the 
ascent.  An  exception  are  bubbles  rising  in  a highly  viscous  medium  (e.g.,  mineral  oil  and 
com  syrup). 

The  results  for  the  spherical  bubbles  only  are  plotted  in  terms  of  the  drag  coefficient 
and  the  Reynolds  number  in  Figure  19,  with  Reynolds  numbers  ranging  up  to  about  400.  The 
drag  curve  for  rigid  spheres  is  also  included.^*  From  it,  the  following  can  be  observed:  The 
drag  curves  of  spherical  bubbles  in  the  various  liquids  fall  between  two  limiting  curves.  As 
upper  limit,  the  drag  curve  of  rigid  spheres  is  obtained,  while  the  lower  limit  is  the  drag  curve 
for  fluid  spheres.  With  decreasing  Reynolds  number,  the  rigid  sphere  curve  connects  with  the 
straight  line  of  Stokes'  Law,  while  the  fluid  sphere  curve  connects  with  the  line  of  Hadamard* 
Rybczynski's  Law.  The  curve  for  the  fluid  spbwes  was  obtained  by  drawing  the  lower  en- 
velope to  the  experimental  curve;  its  accuracy  can  be  confirmed  by  additional  tests  in  other 
liquids  or  by  extension  of  the  theoretical  solution  into  regions  beyond  that  of  very  slow  flow. 

It  will  be  noted  from  Figure  19  that  the  curve  for  mineral  oil,  for  example,  follows  the 
strai^t  line  of  Hadamard-Rybczynski's  Law  over  a certain  region  of  Reynolds  numbers. 

This  indicates  that  the  boundary  conditions  assumed  in  the  analytical  solution  for  fluid 
spheres  are  actually,  fulfllled  and  that  circulation  exists  inside  the  bubble . Circulation  in- 
side bubbles  has  been  observed  experimentally.^^ 

The  experimental  curves  of  Figure  19  also  indicate  an  interesting  aspect  of  the  phe- 
nomenon of  bubble  motion,  namely  that  with  decreasing  Reynolds  number,  the  drag  coefficient 
of  the  bubbles  becomes  equal  to  the  drag  of  rigid  spheres.  This  transition  may  occur  at  a 
Reynolds  number  of  about  40  (as  for  Altered  and  distilled  water)  or  may  not  take  place  until 
very  low  Reynolds  numbers  are  reached,  i.e.,  well  within  the  region  of  slow  flow  (as  for  olive 
oil^^  or  very  viscous  syrup ^^).  Thus,  from  the  experimental  data  available,  it  appears  certain 
that  tiny  air  bubbles  rising  in  any  liquid  follow  Stokes’  Law. 

For  bubbles  behaving  like  rigid  bodies,  thus  indicating  absence  of  motion  inside  the 
bubble,  the  internal  Reynolds  number  (although  nonvanishing)  is  of  no  significance  in  de- 
scribing the  rising  motion  of  the  bubbles.  Likewise,  the  internal  Reynolds  number  cannot  be 
used  to  predict  the  transition  point  at  which  the  drag  of  the  bubbles  becomes  less  than  that 
of  corresponding  rigid  spheres.  Beyond  this  transition  point,  the  internal  Reynolds  number 
might  be  of  importance  in  describing  the  motion  of  the  bubbles. 

Surface  tension  tends  to  make  the  surface  area  of  the  bubble  as  small  as  possible. 

For  a given  volume,  the  conAguration  of  minimum  surface  area  is  a sphere.  This  effect  of 
surface  tension  would  be  most  pronounced  for  bubbles  of  small  radii. 
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ELLIPSOIDAL  BUBBLES 

For  larger  bubble  sizes,  the  surface  forces,  which  are  essential  in  maintaining  the 
spherical  shape  of  a bubble,  become  smaller  in  comparison  to  the  viscous  and  hydrodynamic 
forces,  and  flattening  of  the  bubble  occurs.  This  flattening  to  approximately  an  oblate 
spheroid  results  in  higher  drag  as  compared  to  a sphere  of  the  same  volume.  Figure  20  shows 
the  drag  curves  of  ellipsoidal  and  spherical  cap  bubbles  in  terms  of  the  Reynolds  number. 

The  estimated  extent  of  the  regions  of  ellipsoidal  and  spherical  cap  bubbles  is  indicated 
in  the  figure.  It  will  be  observed  that  tlie  region  of  ellipsoidal  bubbles  for  the  various  liquids 
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Figure  18  - Typical  Shapes  of  Bubbles  of  Several  Volumes  in  the  Various  Liquids 

occurs  at  different  ranges  of  Reynolds  number,  that  for  liquids  of  low  "/I/”  number  a minimum 
in  the  drag  curves  is  reached  at  Reynolds  numbers  of  the  order  of  260,  and  that  ttiese  minima 
occur  near  the  transition  from  spherical  to  ellipsoidal  shape.  Such  minima  are  not  obtained 
for  liquids  of  high  ‘W”  number.  The  drag  coefficients  of  bubbles  in  sqch  liquids  decrease 
until  a constant  value  for  tlie  drag  coefficient  (spherical  caps)  is  attained.  In  the  ellipsoidal 
region  the  curves  are  arranged  according  to  the  magnitude  of  the  “4/”  number,  indicating  that 
the  liquid  properties  contained  in  this  parameter,  namely  surface  tension,  viscosity,  and 
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Figure  20  * Drag  Coefficient  as  a Function  of  Reynolds  Number  for  Ellipsoidal 
and  Spherical  Cap  Bubbles  in  Various  Liquids 


W«b«f  Number 

Figure  21  * Drag  Coefficient  as  a Function  of  Weber  Number  for  Ellipsoidal 
and  Spherical  Cap  Bubbles  in  Various  Liquids 
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density,  are  of  primary  importance  in  the  motion  of  these  bubbles. 

Thus,  in  this  region  the  Reynolds  number  inside  the  bubble  is  of  no  importance  in  the 
description  of  bubble  rise,  since  correlation  was  obtained  in  terms  of  drag  coefficient,  Rey- 
nolds number,  and  **4/**  number. 

In  Figure  21  the  drag  coefficient  of  the  bubbles  was  plotted  as  a function  of  Weber 
number.  It  is  seen  that,  for  the  liquids  tested,  transition  to  a ccmstant  value  of  drag  coeffi- 
cient (spherical  cap  region)  is  reached  at  a Weber  number  of  about  20. 

SPHERICAL  CAP  BUBBLES 

When  the  viscous  and  surface  tension  forces  become  small  relative  to  the  hydrodynamic 
forces,  the  shape  assumed  by  the  bubbles  is  that  of  the  so-called  spherical  caps.  Typical 
shapes  of  these  bubbles  are  shown  in  Figure  18.  The  upper  surface  is  essentially  spherical, 
while  the  lower  surface  varies  from  a highly  irregular  one  for  liquids  of  low  viscosity  to  a 
smooth  surface  for  very  viscous  liquids.  The  configuration  of  the  upper  surface  results  al- 
most exclusively  from  the  hydrodynamic  forces.* 

The  geometric  similarity  of  these  bubbles  was  shown  by  Rosenberg^  and  Davies  and 
Taylor,^ ^ who  determined  a constant  drag  coefficient  of  2.6  for  them.  The  results  of  the 
present  tests  in  a number  of  liquids  confirm  this  value;  Figure  21.  The  velocity  of  spherical 
cap  bubbles  of  given  size  rising  in  any  liquid  can  be  determined  from  the  constant  value  of 
the  drag  coefficient  or  directly  from  the  velocity  curve  (Figure  15).  For  Cq  (8/3)^ f^C/* 

- 2.6,  we  obtain  for  the  rate  of  rise  of  the  spherical  caps  in  all  liquids 

u - i.02y^ 

Thus,  the  velocity  of  rise  of  these  bubbles  is  a function  of  the  bubble  size  only  and 
not  of  the  properties  of  the  liquid  (see  special  case  of  "Dimensional  Analysis"). 

PATH  OF  BUBBLES 

Figures  22-26  show  representative  paths  and  corresponding  shapes  of  bubbles  in  the 
various  liquids.  Three  types  of  motion  of  the  bubbles  were  observed  in  the  experiments: 

(1)  rectilinear  motion,  (2)  motion  in  a helical  path,  and  (8)  rectilinear  motion  with  rocking. 

The  motion  of  spherical  bubbles  is  either  rectilinear  or  helical.  For  ellipsoidal  and  spherical 
cap  bubbles,  all  three  types  of  motion  can  occur.  It  appears  that  the  type  of  motion  may  be 
pre<1icted  from  the  value  of  the  Reynolds  number  at  which  the  motion  takes  place.  Below 
Reynolds  numbers  of  about  800  the  motion  is  rectilinear.  With  increase  in  Reynolds  number 
spiraling  begins  and  increases  in  amplitude  and  firequency  until  a maximum  is  reached.  At 
Reynolds  numbers  of  about  8000,  the  spiraling  disappears  and  only  rectilinear  motion  with 

*Sm,  R*f«rmic* 


Reynolds  Number  59  375  527  594  1110  1960 

Scale:  0.431/1 

Figure  22  - Path  and  Corresponding  Shapes  of  Air  Bubbles  Rising  in  Methyl  Alcohol 
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Figure  23  - Path  and  Conespondiag  Shapes  of  Air  Bubbles  Rising  in  Varsol 


Scale:  0.490/1 

Figure  24  - Path  and  Corresponding  Shapes  of  Air  Bubbles  Rising  in  Cold  Water 
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Figure  25  - Path  and  Corresponding  Shapes  of  Air  Bubbles  Rising  in  Mineral  Oil 
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rocking  is  obtained.  For  the  bubbles  rising  in  mineral  oil  and  the  com  syrup  mixtures,  only 
rectilinear  motion  without  rocking  was  observed.  The  maximum  Reynolds  number  reached 
during  those  tests  was  160. 

The  helical  path  of  the  bubbles  assumes  either  a clockwise  or  counterclockwise  direc* 
tion,  depending  upon  conditions  at  generation.  The  velocity  of  rise  of  these  bubbles  is  not 
affected  by  the  sense  in  which  the  bubble  revolves.  The  major  axis  of  ellipsoidal  bubbles 
is  always  perpendicular  to  the  direction  of  motion. 

The  oscillatory  motion  of  bubbles  is  probably  caused  by  the  periodic  shedding  of  vor- 
tices behind  the  bubble.  Such  vortex  shedding  has  been  observed  experimentally  for  rigid 
spheres  at  the  same  magnitudes  of  Reynolds  numbers  as  for  bubbles.**^ 

BUBBLES  IN  FILTERED  AND  TAP  WATER 

The  results  of  the  experiments  in  filtered  water  at  room  temperature  and  in  cold  filtered 
water  are  presented  in  terms  of  drag  coefficient  and  Reynolds  number  in  Figure  26.  The  re- 
sults of  Bryn,  Allen, ^ and  Gorodetskaya^^  are  also  included.  It  is  seen  that  the  drag  curves 
at  the  two  different  temperatures  coincide  in  the  spherical  and  spherical  cap  region.  In  the 
region  of  ellipsoidal  bubbles,  the  drag  coefficient  at  a given  Reynolds  number  increases  with 
increase  in  “If"  number.  Gorodetskaya  and  Allen  conducted  their  experiments  in  distilled 
water.  A comparison  of  their  experimental  data  with  those  obtained  in  filtered  water  shows, 
within  experimental  accuracy,  no  difference  in  the  drag  of  air  bubbles  rising  in  filtered  and 
distilled  water. 

The  drag  coefficients  for  air  bubbles  rising  in  tap  water  at  two  different  temperatures 
are  given  in  Figure  27.  Gorodetskaya's  results  in  tap  water  at  room  temperature  are  included 
in  the  figure.  For  comparison,  the  drag  curves  for  bubbles  in  filtered  water  at  room  tempera- 
ture and  for  rigid  spheres  are  also  shown.  Again,  in  the  region  of  spherical  and  spherical  cap 
bubbles,  the  drag  curves  at  the  two  temperatures  coincide.  The  value  of  the  minimum  drag 
coefficient  is,  however,  greater  than  that  of  the  corresponding  filtered  water.  In  general,  for 

Reynolds  numbers  up  to  about  800,  the  drag  curves  of  bubbles  in  tap  water  follow  closely  the 

« 

curve  of  rigid  spheres. 

Thus,  the  results  of  the  experiments  (given  in  Figures  5,  26,  and  27)  show  that  for  bub- 
bles (ranging  in  equivalent  radius  from  0.085  to  0.26  cm)  it  is  important  whether  the  motion 
occurs  in  filtered  (distilled)  or  tap  water. 

In  view  of  the  fact  that  merely  filtering  the  water  was  sufficient  to  produce  a change  in 
the  drag  of  the  bubbles,  it  is  indicated  that  the  presence  of  minute  particles  causes  this 
change.  Minute  particles,  most  of  which  can  be  removed  by  filtering,  are  known  to  exist  in 
ordinary  tap  water.  Specifically,  if  such  particles  are  present  in  the  water  a high  concentration 


•It  should  bs  notsd  hors  thst  the  physical  properties  of  tap  water  did  not  differ  from  those  of  fUtered  water. 
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Figure  26  * Drag  Coefficient  as  a Function  of  Reynolds  Number  for  Air 
Bubbles  Rising  at  Their  Terminal  Velocity  in  Filtered  or  Distilled 
Water  at  Different  Temperatures 
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Figure  27  - Drag  Coefficient  as  a Function  of  Reynolds  Number  for  Air  Bubbles 
Rising  at  Their  Terminal  Velocity  in  Tap  Water  at 
Two  Different  Temperatures 
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of  these  psrtioles  would  be  found  at  the  surface  of  the  bubble.  The  particles  at  the  surface 
would  travel  with  the  bubble,  hence  imparting,  in  effect,  a rigid  surface  to  the  bubble.  With 
increasing  bubble  velocities,  the  shear  forces  become  large  in  comparison  to  the  forces  bold* 
ing  the  particles  to  the  surface  and  hence  at  a certain  critical  velocity  no  high  concentration 
of  particles  on  the  surface  can  exist. 

Figure  28  compares  the  drag  curve  of  bubbles  in  tap  water  obtained  in  the  Model  Basin 
study  and  the  experimental  data  of  air  bubbles  rising  in  water  at  room  temperature  as  obtained 
by  various  other  investigators.  Reports  of  these  other  studies  give  no  information  regarding 
the  type  (tap,  filtered,  or  distilled)  of  water  used,  but  presumably  these  experiments  were 
also  conducted  in  tap  water. 

EFFEa  OF  SURFACE-ACTIVE  SUBSTANCES 

The  effect  of  surface-active  substances  on  the  rate  of  rise  of  air  bubbles  has  previous* 
ly  been  investigated  by  Gorodetskaya,^*^  who  added  small  concentrations  of  various  alcohols 
to  water  and  concluded  that,  beyond  a certain  critical  concentration  of  the  surface-active  sub- 
stance, the  rate  of  rise  of  the  air  bubbles  is  not  affected.  Stuke^*^  ran  experiments  with  oxy- 
gen bubbles  rising  in  water  containing  small  concentrations  of  caproic  acid.  The  concentra- 
tions of  the  alcohols  and  the  caproic  acid  were  relatively  small,  hence  the  decrease  in  the 
surface  tension  was  only  about  1 dyne/cm.  In  the  present  study,  the  authors  conducted  tests 
in  water  containing  Glim,  a liquid  detergent.  The  concentration  of  Glim  (0.42  percent  by  vol- 
ume) was  high  enough  to  decrease  the  surface  tension  by  40  dynes/cm.  No  measurable  change 
in  the  viscosity  and  density  of  the  test  liquid  due  to  the  presence  of  Glim  was  noted  (see 
Table  1).  This  was  also  true  for  the  alcohol  and  caproic  acid  solutions.  Results  from  these 
experiments  are  presented  in  terms  of  the  drag  coefficient  and  Reynolds  number  in  Figure  29.* 

The  drag  curve  for  bubbles  in  the  Glim  solution,  as  well  as  the  experimental  data  from 
the  other  investigators  in  water  containing  at  least  the  critical  concentration  of  the  surface- 
active  substance,  follows  the  drag  curve  of  rigid  spheres  to  a Reynolds  number  of  about  200.*  * 
In  the  region  of  Reynolds  numbers  of  10  to  200,  the  drag  curve  for  bubbles  rising  in  a pure 
liquid  having  an  number  very  close  to  that  of  the  Glim  solution***  follows  the  drag 
curve  of  fluid  spheres.  Thus,  the  motion  of  bubbles  in  water  containing  surface-active  ma- 
terials cannot  be  compared  with  that  of  bubbles  in  pure  liquids  on  the  basis  of  drag  coefficient, 
Reynolds  number,  and  “AT*  number,  even  in  the  region  of  ellipsoidal  bubbles.  Although  the 


*Aa  abown  In  th«  Appandix,  no  significant  dlffarance  In  tha  tarminal  valocity  of  oxygan  and  air  bubbles  rising 
in  distillad  wstsr  is  obtainad.  Hanca,  inclualon  of  tha  raaulta  of  tha  taata  with  oxygan  bubblas  in  water  contain- 
ing aurfaca-activa  subatancaa  la  Justified. 

•*For  concentrations  below  tha  critical,  tha  drag  curve  lisa  between  that  of  pure  watar  and  tha  curve  ahown  in 
Plgura 

••♦The  "Af"  nunbar  of  Glim  waa  Z78  X lO'  that  of  Varsol,  for  axanpla,  waa  4.3  X 
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Figure  28  • Drag  Coefficient  as  a Function  of  Reynolds  Number  for  Air  Bubbles 
Rising  at  Their  Terminal  Velocity  in  Water  at  Room  Temperature 
as  Obtained  from  Data  of  Various  Investigators 


Figure  29  • Drag  Coefficient  as  a Function  of  Reynolds  Number  for  Bubbles 
Rising  at  Their  Terminal  Velocity  in  Water  Containing 
Various  Surface-Active  Materials 
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range  of  bubble  sizes  did  not  extend  fully  into  the  region  of  spherical  cap  bubbles,  it  is  quite 
certain  that  the  presence  of  surface* active  substances  will  not  alter  their  rate  of  rise,  which 
was  shown  to  be  independent  of  all  physical  properties  of  the  liquid. 

The  difference  in  behavior  of  the  air  bubbles  must  be  sought  in  the  behavior  at  the 
surface.  A high  concentration  of  molecules  of  the  surface-active  substances  will  be  found 
at  the  surface  of  the  bubble.  As  for  the  case  of  the  tap  water,  these  molecules  would  travel 
with  the  bubble  and  would  impart,  in  effect,  a “rigid*'  surface  to  the  bubble,  that  is  to  say, 
would  impose  the  condition  of  zero  velocity  at  the  boundary. 

Thus,  the  results  of  the  tests  given  in  Figure  29  show  that  these  surface-active  sub- 
stances increase  the  drag  of  the  bubbles  (in  the  region  of  bubbles  having  equivalent  radii  of 
0.08  to  0.30  cm);  beyond  the  critical  concentrations,  any  increase  in  concentration  has  rela- 
tively little  influence  on  the  drag  of  the  bubbles. 

BUBBLES  AS  RldO  BODIES 

In  previous  sections  it  was  shown  that  bubbles  rising  in  pure  liquids*  behave  essen- 
tially like  fluid  bodies  over  a large  range  of  bubble  size,  but  that  below  a certain  critical 
size  (the  size  being  different  for  various  liquids),  the  bubbles  behave  like  rigid  bodies,  that 
is  to  say,  the  drag  of  the  bubbles  equals  that  of  corresponding  rigid  bodies. 

A possible  explanation  of  the  anomaly  of  behavior  of  the  gas  bubbles  is  given  by 
Boussinesq’s  dynamic  surface  tension.^  As  pointed  out  under  “Theoretical  Solutions,*’ 
for  small  bubbles  the  effect  of  the  dynamic  increment  of  surface  tension  increases  the  drag 
to  the  value  of  corresponding  rigid  bodies.  With  increase  in  bubble  size,  this  effect  becomes 
negligible  and  the  drag  of  the  bubble  equals  that  of  a fluid  body.  The  dynamic  increment  in- 
cludes a constant  of  proportionality  (surface  viscosity)  which  is  a function  of  the  two  fluids 
composing  the  interface.  Therefore,  for  air  bubbles,  the  transition  region  from  “rigid**  to 
fluid  bodies  would  be  different  for  the  various  liquids.  There  is,  however,  no  experimental 
evidence  that  dynamic  surface  tension,  as  postulated  by  Boussinesq,  exists. 

From  a “hydrodynamic**  point  of  view,  the  reason  for  the  transition  of  the  bubbles  to 
“rigid**  bodies  is  not  clear.  As  mentioned  earlier,  the  mere  inclusion  of  surface  tension  as 
a pressure  drop  in  the  boundary  conditions  does  not  alter  the  analytical  solution  for  fluid 
spheres.  Hence,  it  appears  that  the  presence  of  surface  tension  should  have  no  effect  on  the 
motion  of  the  bubble,  except  in  maintaining  the  spherical  shape.  Thus  the  anomalous  behavior 
of  the  bubbles  cannot  be  explained  in  terms  of  “hydrodynamics,**  but  must  be  sought  in  terms 
of  a surface  phenomenon.  If  it  could  be  shown  analytically  (if  only  in  the  region  of  slow  flow) 
that  equality  of  drag  of  corresponding  rigid  bodies  and  bubbles  also  implies  equality  of  bound- 
ary conditions  at  the  surface,  then,  as  in  the  case  of  rigid  bodies,  the  velocity  of  the  entire 

*Mixtur«s  Buch  th«  13  pcrcant  Bthyl  aleohol-w«t«r  mlstura  ar*  included  in  thia  categoiy. 
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flivfaoe  of  the  bubble  must  vanish.  The  surface  must  then  be  able  to  hold  molecules  of  the 
pure  liquid,  just  as  in  the  case  of  tap  water  and  surface-active  substances  the  surface  aU 
tracts  and  holds  a high  concentration  of  particles  or  molecules  of  the  surface-active  sub- 
stances. The  molecules  at  the  surface  would  travel  with  the  bubble  and  hence  would,  in  ef- 
fect, give  the  same  boundary  conditions  as  a rigid  surface.  As  the  shear  forces  become 
larger  in  comparison  to  the  forces  holding  the  molecules  at  the  surface,  "rigidity"  at  the 
surface  cannot  be  maintained;  circulation  inside  the  bubble  ensues  and  the  drag  of  the  bubble 
becomes  smaller  as  compared  to  that  of  a rigid  body. 

In  addition,  since  it  was  not  possible  to  correlate  the  results  of  the  experiments  on 
the  motion  of  bubbles  in  the  gravity  pressure  Held  in  terms  of  nondimensional  parameters 
formed  from  the  usual  liquid  properties  (viscosity,  surface  tension,  density),  further  work  on 
freely  rising  bubbles  is  necessary  before  the  results  obtained  from  such  tests  can  be  utilized 
and  before  the  more  complicated  behavior  of  bubbles  in  variable  pressure  gradients  can  be 
understood.  Particularly,  an  understanding  of  the  reason  for  the  transition  of  bubbles  from 
fluid  to  "rigid"  bodies  as  well  as  a criterion  for  this  transition  is  most  desirable,  since  such 
transition  might  be  influenced  by  the  magnitude  of  the  pressure  gradient.  In  the  region  of 
Reynolds  numbers  where  the  bubbles  behave  like  rigid  spheres,  the  pressure  gradient  probably 
has  no  effect  on  the  drag  coefficient. 


SUMMARY 

As  the  size  of  the  bubbles  was  increased  in  the  tests,  a change  in  bubble  shape  from 
spherical  to  ellipsoidal  to  spherical  cap  shape  was  observed  in  all  liquids.  The  volumes  at 
which  these  transitions  occur,  however,  varied  with  the  properties  of  the  liquid.  For  spheri- 
cal bubbles  of  given  volume  the  results  show  that  the  viscosity  of  the  liquid  is  the  most  im- 
portant property  determining  the  rate  of  rise.  For  ellipsoidal  bubbles,  the  surface  tension 
assumes  greater  importance.  Spherical  cap  bubbles  rise  independently  of  the  properties  of 
the  liquid. 

The  results  show  that  the  motion  of  air  bubbles  rising  at  their  terminal  velocity  in  a 
gravity  field  cannot  be  described  completely  by  use  of  dimensionless  parameters  formed  from 
the  usual  liquid  properties  (viscosity,  surface  tension,  density),  the  equivalent  radius  of  the 
bubble,  the  acceleration  of  gravity,  and  the  terminal  velocity. 

The  drag  coefficients  of  tiny  spherical  bubbles  coincide  with  those  of  corresponding 
rigid  spheres.  With  increase  in  bubble  size,  a decrease  in  the  drag  as  compared  to  that  of 
rigid  spheres  occurs  in  some  liquids.  This  change  in  the  drag  is  due  to  the  development  of 
.circulation  inside  the  bubble.  The  drag  curves  of  the  spherical  bubbles  rising  in  various 
liquids  fall  between  two  limiting  curves,  namely  the  drag  curves  of  rigid  and  fluid  spheres, 
respectively.  It  was  not  possible  to  determine  a criterion  for  the  transition  region  of  the 
bubbles  from  "rigid"  to  fluid  spheres. 

The  region  of  ellipsoidal  bubbles  extends  over  different  ranges  of  Reynolds  numbers 
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for  the  various  liquids.  For  liquids  of  low  number  (say  less  than  10**^),  a minimum  in 
the  drag  curve  is  reached  at  Reynolds  numbers  of  the  order  of  260.  These  minima  occur  near 
the  transition  from  spherical  to  ellipsoidal  shape.  Such  minima  in  the  drag  curve  are  not  ob- 
tained for  liquids  of  high  “/!/"  number.  For  the  liquids  used,  transition  to  spherical  caps  is 
completed  at  a Weber  number  of  about  20. 

The  drag  coefficients  of  spherical  cap  bubbles  ate  independent  of  bubble  size  and 
have  a constant  value  of  2.6-  The  rate  of  rise  of  these  bubbles  as  a function  of  the  equivalent 
radius  is  given  by  the  experimentally  determined  relation: 

V - 1.02 

For  bubbles  (ranging  in  equivalent  radius  from  0.08  to  0.26  cm)  rising  in  tap  water,  an 
increased  drag  as  compared  to  bubbles  in  clean  (filtered  or  distilled)  water  was  observed. 

The  presence  of  certain  surface-active  substances  in  the  water  similarly  increases  the  drag 
of  bubbles  (ranging  in  equivalent  radius  from  0.08  to  0.80  cm)  as  compared  to  bubbles  in 
pure  water.  Beyond  a certain  critical  concentration  of  these  surface-active  substances,  an 
increase  in  concentration  has  relatively  little  influence  on  the  drag  of  the  bubbles. 

Tests  to  determine  the  effect  of  the  container  walls  on  the  velocity  of  rise  indicate 
the  absence  of  such  effect  for  the  range  of  bubble  volumes  and  container  sizes  tested. 
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APPENDIX 

RATE  OF  RISE  OF  GAS  BUBBLES  IN  DISTILLED  WATER 

It  is  of  interest  to  compare  the  rate  of  rise  of  bubbles  of  different  gases  in  distilled 
water.  As  pointed  out  under  "Previous  Experimental  Work,*’  a number  of  tests  have  been 
conducted  using  oxygen  bubbles.  In  Figure  80,  the  velocity  of  sir  and  of  oxygen  bubbles  in 
distilled  water  is  given  as  a function  of  the  equivalent  radius  of  the  bubble.  The  results 
indicate  no  signiScant  difference  in  the  rate  of  rise  of  the  bubbles  composed  of  oxygen  and 
those  of  air. 


Figure  80  • Terminal  Velocity  of  Gas  Bubbles  in  Distilled  Water 
as  a Function  of  Bubble  Sixe 
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